Spin alignment of vector mesons in heavy ion and proton - proton collisions 
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The spin alignment matrix element poo for the vector mesons K*^ and 0(1020) has been measured 
in RHIC at central rapidities. These measurements are consistent with the absence of polarization 
with respect to the reaction plane in mid-central Au + Au collisions whereas, when measured with 
respect to the production plane in the same reactions and in p + p collisions, a non-vanishing and 
px-dependent poo is found. We show that this behavior can be understood in a simple model of 
vector meson production where the spin of their constituent quarks is oriented during hadronization 
as the result of Thomas precession. 

PACS numbers: 25.75.-q, 13.88.+e, 12.38.Mh, 25.75.Nq 



The study of spin polarization of produced hadrons in 
reactions at high energies has opened a window to the un- 
derstanding of the underlying dynamics of quark recom- 
bination. In the context of heavy-ion collisions, polariza- 
tion studies can also help to understand the evolution of 
the system from its early stages [H, 0: ■ 

Polarization analyses require to determine a given di- 
rection that serves as the spin quantization axis. From 
the experimental point of view, it is possible to determine 
two directions: the normal to the reaction and the nor- 
mal to the production planes. The fist plane is defined as 
the one containing the impact parameter and the beam 
direction vectors whereas the second one is defined as 
containing the hadron's final momentum and the beam 
direction vectors. 

Recently, the STAR collaboration has reported mea- 
surements of the spin polarization properties of the vec- 
tor mesons (f> and K* [J]. These measurements refer to 
the 00 component of the so called spin alignment density 
matrix p, which is the density matrix for a two-spin one- 
half system in a triplet state, expressed in terms of the 
coupled basis 0]. Recall that a value p°" = 1/3, means 
that the spin of the vector meson is not aligned with re- 
spect to the chosen quantization axis. Deviations from 
this value indicate a degree of polarization of the vector 
spin which ultimately might reflect a polarization of the 
constituent quarks. 

The experimental findings reported can be summarized 
as follows: When the spin alignment is referred to the 
reaction plane in Au + Au collisions at ^/s'nn = 200 GeV 
and measured at mid-rapidity, poo — 1/3, and it remains 
constant both as a function of p± in the range < p± < 5 
GeV, for mid-central collisions, and as a function of the 



average number of participants in the same p± range, for 
both (f) and K*. When the spin alignment is referred 
to the production plane and measured at mid-rapidity, 
both for p -|- p and mid-central Au -|- Au collisions at 
y/SNN ~ 200 GeV, Poo > 1/3 and it has a concave shape 
as a function of p± in the range < p_L < 5 GeV with 
minima at slightly different intermediate values of p± for 
and K*. 

An interesting observation that can be inferred from 
the above listed results is that the dynamics of hadron 
formation seems to play a role to orient the spin of quarks 
that form meson. 

A possible origin of a quark polarization driving the po- 
larization of vector mesons is discussed in Refs. [l| which 
study the transfer of local relative angular momentum in 
peripheral nuclear collisions to the quark spin polariza- 
tion by means of rescattering during the reaction. How- 
ever, this mechanism also predicts a small global polar- 
ization growing almost linearly with impact parameter 
which, according to the aforementioned results is not ob- 
served in data. 

Another interesting possibility that has only been stud- 
ied in the context of hyperon polarization 0] is the sce- 
nario where the spin of a quark is oriented during the re- 
combination process. The semiclassical picture account- 
ing for the polarization is the Thomas precession pro- 
duced by the accelerating force that pulls a slow moving 
quark {q'^ ) to form a fast moving hadron Q . This mech- 
anism also predicts that if the quark is fast (g-^) and is 
decelerated to form the hadron, its polarization will be 
of opposite sign compared to the case when it is acceler- 
ated. The pulling force is required to not be parallel to 
the original quark velocity since the Thomas frequency is 
a vector formed by the cross product of the force F and 
this velocity /3, namely 
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where 7 is the Lorentz gamma-factor. The polarization 
is given by Q 



UJ. 



AE' 
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u)!^'^ is the 



where the — (+) sign refers to the q" (q^) 
magnitude of the Thomas precession frequency for and 
q-^ , respectively and AE is the change of energy in the 
process of hadron formation. 

In this work we use the Thomas spin precession mech- 
anism to describe the spin alignment of vector mesons 
produced at central rapidity in Au -I- Au and p 4- p col- 
lisions at y^SNN = 200. We show that under very simple 
assumptions, data for pqq are well reproduced within this 
approach. 

The physical picture we use is that of a fast quark 
that decelerates and a slow one that accelerates to form 
a fast moving hadron. In the process, Thomas precession 
makes the spin of the former to acquire a positive polar- 
ization whereas the latter acquires a negative one. For 
central rapidities, the large momentum component of the 
hadron will thus be its transverse momentum whereas 
the small component will be its longitudinal one. . We 
will assume that in the beam collision, a hard interaction 
produces a fast quark moving with a large transverse mo- 
mentum and, to simplify matters, a vanishing longi- 
tudinal momentum. This fast quark combines with the 
slow one, that we assume moves originally mainly in the 
longitudinal direction with momentum pjj and, also for 
simplicity, take it with vanishing transverse momentum. 
This sharp difference in the original direction of motion 
of the recombining quarks is at the core of the produced 
polarization since, as we proceed to show, it gives rise to 
a distinct p± dependence of poo which seems to be also 
observed in data. 

In order to form the hadron, which should move with 
an intermediate value of momentum, between that of 
the q-f and of the q**, the fast quark should slow down 
whereas the slow quark should speed up. Notice that 
for this mechanism to work, there is no need to assume 
that the process happens only in cither a proton-proton 
or a nucleus-nucleus collision. Notice also that the mo- 
mentum of the formed hadron provides a fixed direction 
to define that q" (g-^) decelerates (accelerates), whereas, 
when referred to the reaction plane, no such fixed direc- 
tion exists, since the direction of the impact parameter 
vector changes from one reaction to another and in such 
situation either quark can accelerate or decelerate. 

The pulling force is equal to the change in momentum 
Ap of the given quark, in the interval of time At for the 
recombination process to happen, that is 



average over this time interval [a| , namely 
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where in the last line we have changed the time average of 
the sine of the angle by the average sine of the angle be- 
tween the quark velocity vectors and their corresponding 
change in momentum. From Eqs. ([2]) and we see that 
the calculation of the q^ and q'^ polarizations reduces to 
computing the magnitudes of their change in momentum 
Ap'''^ , their (sin 6')'''''' and the change in energy AE. 

We first compute the change in momenta. For q" , we 
have 



Ap^ = y'(p;;/^-pp2 + (pVif)2^ 



(5) 



where denote the parallel (transverse) component 

of the g** in the hadron H. Let us assume that, in order to 
enhance the recombination probability, the rapidity of g^ 
has to be within the hadron's one. Under this assumption 
we can write 
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where we have set TOj 
can write 



m since 



Pi 



0. Therefore, we 
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where we have introduced the definitions for the momen- 
tum fractions that the g^ has inside the hadron. 
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and have neglected the longitudinal hadron's momentum 
with respect to its transverse one. 

Similarly, for the q^ change of momentum when decel- 
erating, notice that we have 



Jlii 



J/H 
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Ap 
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Thus LOT for the given quark can be computed as the 



Since (3^ initially points along the perpendicular direc- 
tion and, although this velocity vector changes so that 
the final hadron's momentum eventually picks up a lon- 
gitudinal component, the dominant component of the qf 
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in the hadron is the transverse one. Therefore, for the 
cross product of Ap^ with (3-^ we get 



Ap^ X (3^' = pi' 13^ (si 



= (1 -.T||)pf (sin^ 
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where we have enforced momentum conservation xj_ + 
a;|| = 1 and have approximated ~ 1. 

Now, we proceed to compute the average sine of the 
angles. To compute (sin0)*, notice that given that ini- 
tiaUy /3* points afong the parallel direction and Ap" is 
basically directed along the perpendicular direction, the 
initial angle between these vectors is 7r/2 and the aver- 
age one should be close to 7r/4. Thus (sin 6')" ~ l/\/2. 
Rather than approximating [3'^ (and therefore 7^*), we in- 
troduce a factor a for this polarization and let this vary 
such that < a < 1. Thus we write 
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To compute (sin 6) ^ , notice that since Ap-^ and (3^ are 
almost perpendicular, (sin^)-^ ~ 1. 

The change in energy is common to both the acceler- 
ating and the decelerating qf 

AE = {[{p{f + {plf + (mff]'/' 

+ [{plf + iplf + im^V 

C {[(p{)2 + (™/)2]l/2. 

- iip^)' + im^n/'} 

[mff 
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where we have set = = 0, neglected [p^Y com- 
pared with (px)^ ^^cl expanded the square roots as- 
suming that the transverse momentum components are 
large. Introducing the relation between the hadron and 
q^ transverse momenta = p^j with < z < 1, we 
get 
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where we have made use of the assumption that the ra- 
pidity of the coincides with that of the hadron. Notice 
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FIG. 1: (Color online) poo as a function of p± for (^(1020) 
using the model parameters described in the text, compared 
to data from STAR for p -|-p and Au -|- Au collisions at 
centrality 20-60 %, measured with respect to the production 
plane. For clarity, the p± for p -I- p data has been displaced 
by 0.09 GeV with respect to the reported central value. The 
statistical and systematic errors have been added in quadra- 
ture. For comparison, we also draw the constant value 1/3 
that represents the absence of polarization. 



that the above expression can be simplified. In particu- 
lar, since p^ = s'mhy^ , we can write 
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(12) Therefore. AE can be expressed as 
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We enphasize that the approximation to compute AE is 
such that Eq. ([TS]) is valid for p± > and that for 
z < 1 the validity of the approximation can be extended 
to lower values of p± . 

It is now easy to compute the polarization for the slow 
and fast quarks by means of Eq. ([2]) and from them, the 
Poo density matrix element given by 



Poo — 



(16) 



Let us first proceed to apply the model to compute poo 
for (f>. This case is the simplest one to treat within our 
approach since the quark content of this particle is ss and 
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thus cither one of these quarks can be thought of as being 
the fast (g^) or the slow {q'^) one. Rather than making 
an exhaustive search in the parameter space, we choose 
reasonable values for them. We first fix the 4> and strange 
quark masses to be m"^ = 1.02 GeV, m^'^ = 0.5 GeV. 
The rapidity value we use is = 1 and the formation 
time At = 1 fm. For the fractions of longitudinal and 
transverse momenta that the slow quark has inside the 
(j) we take xj^ — ~ 0.5. The fraction of the transverse 
momentum carried by the from the fast quark is taken 
as z = 0.9. 

Figure [1] shows Ppg as a function of compared to 
data from STAR [J| for p + p and Au + Au collisions 
at centrality 20-60 %. A good description is obtained for 
a = 0.25. 

We now proceed to apply this model to the case of 
K* , whose quark content is ds. However, in this case 
one needs to be careful since the symmetry between the 
masses, present in the description of (f>, is absent. Conse- 
quently the spin alignment has to be treated in average. 
To this end, let us first take a simple scenario and con- 
sider the arithmetic average in the following manner 



Poo " 2 ^^00 



/=s, s=d 



f=d. s=s 
Poo 
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Figure [2] shows Poh as a function of compared 
to data from STAR [i| for p + p and Au + Au colli- 
sions at centrality 20-60%, measured in the production 
plane. The curves are computed such that we employ the 
same set of parameters as in the computation of (with 

nA = 0.3 GeV) except for the value of z. The reason is 
that, whereas in the case that q-^ = s one can think that 
in order for this fast quark to pick up a d, the momenta of 
s and K* are similar, when q^ = d, its momentum must 
be larger, given the mass difference between K* and d. 
Thus if = s we choose z = 0.9 whereas when q-l" = d 
we use z = 0.3. The upper dashed curve represents the 
=s, s=d 



case for Pqq * whereas the lower dashed curve is for 

Poo'^ The intermediate solid curve represents Pqq 

as the algebraic average of the above, as in Eq. (fT7|) . 

An alternative approach is to consider that the Pqq 
can be computed by the substitution 



=d 



(18) 



that is, by the average product of polarizations, in 
Eq. (fT6|l . Figure [2] shows also this possibility represented 
by the intermediate dotted curve, using the same set of 
parameters for the cases were q-^ = s, q'' ^ d and q^ ~ d, 

= s, as discussed above. As can be seen from the fig- 
ure, no significant difference is found with the case where 
the average is taken with the functions poo and both ap- 
proaches give a good description of data. 

In conclusion, we have shown that data on the vector 
mesons (j) and K* spin alignment with respect to the pro- 
duction plane in Au 4- Au and p -I- p collisions are well 
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FIG. 2: (Color online) poo as a function of p± for K* using 
the model parameters described in the text, compared to data 
from STAR for p -|-p and Au -I- Au collisions at centrality 
20-60 %, measured with respect to the production plane. For 
clarity, the p± for p + p data has been displaced by 0.09 GeV 
with respect to the reported central value. The statistical and 
systematic errors have been added in quadrature. The upper 

dashed curve represents the case for PqcT^' whereas the 
lower dashed curve represents the case for Pqq'^' The in- 
termediate solid curve represents as the algebraic average 
of the above. The intermediate dotted curve represents the 
case where poo is computed using the average product of po- 
larizations. For comparison, we also draw the constant value 
1/3 that represents the absence of polarization. 



described by assuming that these hadrons are produced 
by the recombination of a slow and a fast quark that in 
the process become polarized in opposite directions due 
to Thomas precession. In this plane, the momentum of 
the hadron provides a fixed direction to define whether 
a valence quark accelerates or decelerates. The mecha- 
nisms also clarifies the fact that when the spin alignment 
is referred to the reaction plane, poo vanishes, given that 
no such fixed direction exists, since impact parameter 
vector changes from one reaction to another. 

In the near future ALICE at the LHC will have the 
capability to measure and reconstruct (f> and K* mesons 
with larger statistics Q. In addition, its particle identi- 
fication will allow these meson's p± to be measured be- 
yond 5 GeV, well into the region where energy losses 
become important and also where fragmentation (as op- 
posed to the recombination picture we are using here), 
becomes the dominant particle production mechanism. 
It will thus be interesting to study how the polarization 
of vector mesons changes when including these effects. 
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This is work for the future. 
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